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We extend previous results in spherically symmetric gravitational systems coupled with a massless 
scalar field within the loop quantum gravity framework. As starting point, we take the Schwarzschild 
spacetime. The results presented here rely on the uniform discretization method. We are able to 
minimize the associated discrete master constraint using a variational method. The trial state for the 
vacuum consists of a direct product of a Fock vacuum for the matter part and a Gaussian centered 
around the classical Schwarzschild solution. This paper follows the line of research presented by 
Gambini, Pullin and Rastgoo [T] and a comparison between their result and the one given in this 
work is made. 
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Introduction 

Loop quantum gravity (LQG) is a promising proposal for the canonical quantization of general relativity [2]. It 
provides a non-perturbative and mathematically rigorous description of the kinematical sector. Nevertheless, as it is 
well known, the treatment of the dynamics of the theory remains an open problem. The main difficulty of this task 
is due to the fact that the constraints in general relativity do not satisfy a Lie algebra. 

During the last years, there have been advances in the understanding of the spherically symmetric gravity within 
the LQG framework, both in the vacuum and in the case coupled with matter [lj [3] . Gambini, Pullin and Rastgoo 
studied the problem of a gravitational system with spherical symmetry coupled to a spherically symmetric massless 



2 



scalar field in a discretized space. They employed the uniform discretization method [IHHj to deal with the fact that 
the Hamiltonian constraint satisfies a non-Lie Poisson bracket with itself (after a gauge fixing of the diffeomorphism 
constraint). This uniform discretizations technique provides a method to quantize totally constrained systems. It is 
as the 'master constraint' program [7J an approach to find an alternative for the Dirac quantization used in LQG. 

More specifically, in [T] the scalar field starting from the Minkowski spacetime is studied. The authors performed 
a minimization of the master constraint through a suitable variational method. In this way, one can check whether 
the theory provides us with the vacuum state of the scalar field centered around a Minkowskian background. In the 
present paper we extend that work to the Schwarzschild case. Although it is a modest objective, our motivation is 
twofold. On the one hand, we want to test the methodology presented by Gambini et al. in a less trivial case and, 
on the other hand, we want to compare the final result for the Schwarzschild spacetime with the one in Minkowski. 
As we are going to show in the paper, although in this case the procedure becomes more involved than the previous 
one, the final result is, at the leading order, independent of the Schwarzschild radius. 

This paper is organized as follows. In section [T] we review the uniform discretizations method, that will be the 
arena for our study. In section [Tl] we summarize the classical treatment for spherical symmetric gravity coupled 
with a massless scalar field. More concretely, we will write explicitly the Hamiltonian constraint for this system. In 
the next section III we construct the trial state (composed by a gravitational part quantized polymerically and the 
Fock vacuum for the matter part) that we will use later in section IV to evaluate and minimize the discrete master 
constraint. In order to make the main text and ideas more clear, we have gathered most of the technical calculations 
in the appendix and we added the corresponding references where necessary. Of special importance is the appendix 
|B} where we comment the solutions of the Sturm-Liouville equation that plays a central role for the Fock quantization 
of the matter part. 



I. UNIFORM DISCRETIZATIONS 



The uniform discretizations technique [1H6] is a method to quantize totally constrained systems. It is, as the master 
constraint program [7J, an approach to find an alternative for the Dirac quantization used in canonical LQG. Although 
LQG showed remarkable results, especially concerning microscopic pictures of the early universe and black holes, the 
dynamics of the full theory still faces problems. One of those is that the quantum constraint algebra fails to portray 
the classical algebra exactly. Among the other options under research is the uniform discretization method. 

The quantization of discretized, totally constrained systems faces the problem that, in general, the constraint 
algebra fails to close, i.e., the systems turn out to be second class even if the continuum theory is first class. If one 
now uses the Dirac method to quantize the discrete theory, the degrees of freedom will not match the ones of the 
continuum. Furthermore, in general the symmetries of the theory are not recovered and taking the continuum limit 
one will not find the quantized version of the theory one started with. 

If one has first class systems the uniform discretization technique is akin to the Dirac procedure, so we concentrate 
on second class systems here. 

Starting with the study of a classical, totally constrained system with N configuration variables and M constraints 
<f>j\ we consider the discretization of the constraints with a parameter e, such that <f>j = lim e _>.o </>j • As mentioned 
beforehand we assume the constraints to be second class in the discrete theory: 

{^Al} = q k m ^ m + A%. (i) 

In order to get first class constraints in the continuum theory we impose lim c ^.o Aj k = , and for the structure 
functions of the continuum theory C jfc m = lim e ^o Cj k m . 
We now construct a discrete "master constraint" 

i M 

h^Je^) • (2) 
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We proceed by introducing the time evolution for a discrete variable A n as 

A n+1 = eS ■ M 'U n := A n + {A n , ff} + \ {{A n , H e } , ff} + . . . (3) 

It is obvious that H e is a constant of motion under this time evolution. We fix its value to H c = S 2 /2 and define the 
quantities Af :— 4>\/5. With these definitions the evolution of the constraints is given by 



$(n + 1) = $(n) + q k m <f m AlS + A%X%6 + 0(6 2 ) 



(4) 
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Taking the limits <5 — >• and £->Owe recover the evolution equations for the constraints in the continuum theory 

6 e An + 1) - 6 e An) 

& = Km ^ > 6 VA = C jk m cf> m X k , (5) 



where A^ = lim^o Af are Lagrange multipliers. Let us remark that although at first sight the A^ seem not to be free 
in the discrete theory, they are only determined by the constraints evaluated on the initial data, and thus can be 
chosen arbitrarily. 

Applying the Dirac procedure to the discrete theory would lead us to different degrees of freedom in the continuum 
theory, due to the fact that the constraints in this case are second class. This can be shown in the following way. 
Within the Dirac procedure one would impose the constraints 0| and construct the total Hamiltonian Ht = 
with a number M of Lagrange multipliers Cj. Consistency together with the discrete constraint algebra ([I]) would 
require 

&A< jk = 0. (6) 

Therefore the Cj would generally not be free but restricted. This implies that one would have more observables than 
in the continuum theory, i.e., more than 2N—2M, so the degrees of freedom of the discretized theory would not match 
the ones of the continuum theory. Thus, the uniform discretization method is better suited for studying a discretized, 
totally constrained system. 

The uniform discretization approach provides a method to handle the fact that the constraints become second class 
through discretization, without changing the number of degrees of freedom. Moreover the discrete constraint algebra 
mimics the one of the continuum. 



II. SPHERICALLY SYMMETRIC GRAVITY COUPLED TO A SCALAR FIELD 



In this section we review the treatment of the spherically symmetric gravity in Ashtekar variables. This symmetry 
reduction was first performed by Bengtsson [8]. Later on, the definition of the invariant connection components was 
slightly changed to ensure correct transformation properties [S]. We will use the invariant connection formulation 
based on the deep investigation on symmetry reduction of connections in [10] . Using the invariant formulation of the 
basic variables we find the constraints in the spherical symmetric framework. Furthermore we will couple a spherical 
symmetric massless scalar field to gravity, following pQ. Subsequently we perform a gauge fixation in order to solve 
the diffeomorphism constraint. Since the spherical symmetric theory is manifestly gauge invariant we are left with 
the Hamiltonian constraint. The variational method we are using is based on the uniform discretization method and 
thus we end this section with the discretization of the Hamilton constraint. 

In order to describe spherically symmetric spacetimes we use spherical coordinates (r, ip, 9) on a spatial manifold 
with topology £ = M + x S 2 . The invariant connection, densitized triad and extrinsic curvature can be written as: 

A = A r A 3 dr + (AiAi + A 2 A 2 ) dO + ((AiA 2 -A 2 Ai) sin 9 + A 3 cos 9) dtp, (7) 
E = E r sin 6A 3 d r + (E 1 A 1 + E 2 A 2 ) sin 6 d e + {E 1 A 2 - E 2 A l ) d v , (8) 
K = K r A 3 dr + {K 1 A 1 +K 2 A 2 ) d9 + ( KiA 2 -K 2 A 1 ) sin 9 dip, (9) 

where A;, E l ,Ki (i = r, 1, 2) are arbitrary functions on K + and A] are generators 1 of su(2) [TUl 111) . 

In order to get canonically conjugate variables and suitable fall off conditions we use the extended ADM phase 
space variables [HI [12] . For convenience we introduce the following change of variables 

K v {r) = V 'K^r? + K 2 (r) 2 , E v (r) = \J i? 1 (r) 2 + E 2 (r) 2 , (10) 
that implies the Poisson relations 

{K r (x),E r (x')} =2G5(x-x'), (11) 

{K v (x),E«>{x')} = G6(x-x'). (12) 



Aj = 



—icrj/2, where oj denotes a rigid rotation of the Pauli matrices. 
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In terms of the canonical conjugated variables {K r , E r , K v , E v } we can write the constraints as: 

g = o, 

GC = K'E*-\K r (E r )', 



GH = 



'E r {E r )" VE r (E r y(E'*y T ^ t {(E r )') 2 E*K 2 E* 



2Ev 



2{E*) 2 



- y/E r K r K v -r : — — 



(13) 
(14) 

(15) 



where Q, C and H stands for the Gaufi, diffcomorphism and Hamiltonian constraint respectively; G is the Newton 
constant; and /' denotes the derivative of the function / with respect to the radial coordinate. Notice that the 
variables are manifestly gauge invariant. 

Now we add the constraints corresponding to the matter content (the massless scalar field <f>) [15] : 



2EVy/\E F 



((P r ) 2 (0') 2 + (-P*) 2 ) 



(16) 
(17) 



Then, the Hamiltonian and diffcomorphism constraint of the reduced theory minimally coupled to the massless 
scalar field are 



H = 



1 

G 



Mg E*Kl (|gr|/)2 ^\W\\E^\'(E^y VW\(\E r \)" 
2^\W\ ZVW\ 8VW\E<p 2(^) 2 r 2Ev 

(P^) 2 \E r \ 3 ^ 2 ((t) r ) 2 



2^/\W\E^ 2E<p 



(18) 
(19) 

(x + a) 2 as in [3 , where a; is a radial coordinate and a is a constant 
that we will relate with the Schwarzschild radius later on. Then we solve the diffeomorphism constraint, which leads 
to 



C = UE*K' v -\{E r )'K r 



Now, we proceed by gauge fixing E r 



r r 

Rescaling the lapse function N — > NG2r/E v one gets for the Hamiltonian constraint 

H = H v 



vac ~T" ( ' I , 

-a)K\ 

P'l {x + a)\fi) 2 (x + a)K v Pt<t>' 



-x-(x + a)K 2 + 

52 



2 (x + a) 3 \' 
' (E<p) 2 



(E*) 2 



{Ev) 2 



Ef 



(20) 

(21) 
(22) 

(23) 



Now, we proceed to discretize and polymerize the Hamiltonian constraint. We consider 



r — » r(i) , H - 



(r) -> 



p<) 



P 



+ 1) - m 



For the sake of simplicity, we will use a constant parameter p for the polymerization, like in the early loop quantum 
cosmology (LQC). 

Finally, the expression for the discretized Hamiltonian constraint is 



H(i) =- (1 - 2A)e + r(i + 1) 



sbi 2 (pKji + l)) 



r(i) 



sm 2 (pK v (i)) r(i + l) 3 e 2 r{ife 2 



Ev(i + l) 2 Ef{i) 2 



r(i) 4 



e £f +^(^+D-^(or-2 



-j^r^ m% + 1) - H*)) - epv 



(24) 
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where we are employing units such that G = £ 2 , c = h — 1. We have also introduced a cosmological constant A which 
is related to the vacuum energy of the scalar field p vac = |tA. We need to introduce this term in order to consider 
the Fock vacuum of the scalar field. This vacuum energy term comes from the fact that we impose the expectation 
value of the matter constraint in the semiclassical limit £ p — > to vanish. The value of p V ac can then be evaluated by 
introducing the Fock vacuum in a fixed background setting (we will choose our vacuum to be the one defined in the 
asymptotic region of the Schwarzschild spacetime) and computing the "energy of the vacuum", i.e., the expectation 
value of the matter Hamiltonian. We find that the vacuum energy at the order of interest is 

7T 

Pvac = ' 

As pointed out in [T], the introduction of the cosmological constant has some disadvantages. Introducing a cosmo- 
logical constant before the symmetry reduction implies ending up with a non constant term in two dimensions, which 
cannot give rise to the vacuum energy. Because of this, the theory does not stem from a dimensional reduction of the 
full theory if one introduces the cosmological constant in the way we did. To circumvent this problem one can use the 
fact that, unlike in the full four dimensional theory, there is already a constant term in the Hamiltonian constraint. 
The cosmological constant at that level can be interpreted as a rescaling of the radial coordinate. However this 
rescaling has the disadvantage that the volume of spheres is not 4irR 2 anymore, i.e., the full theory one approximates 
has topological defects. 



III. CONSTRUCTION OF QUANTUM TRIAL STATES 



Although the discrete Hamiltonian constraint (24 1 fails to close a first class algebra, it can be shown that using 
the uniform discretization approach one can consistently treat the problem by minimizing the resulting Hamiltonian 
constraint. In order to achieve that, we use the variational technique described in [TJ. 

We have the opportunity to consider the gravitational sector and the matter sector separately, since classically the 
scalar field vanishes in the vacuum. We thus construct a polymeric Hilbert space % V o\y for the gravitational variables. 
The gravitational trial state \^ff) is then constructed as a Gaussian in phase space with width a centered around the 
Schwarzschild solution. With the help of this trial state we find an effective Hamiltonian for the matter variables. Out 
of the equations of motion for this Hamiltonian we construct solutions based on creation and annihilation operators. 
Thus we introduce the Fock vacuum state |0) for the matter sector. In this way, we make contact with the usual 
treatments of quantum field theory in curved spacetimes. 

The trial state used for the minimization of the master constraint is then constructed as a direct product of the 
gravitational Gaussian state and the Fock vacuum 

|^ rial ) = |**) ® |0) . (25) 



A. Gravitational part of the trial states 



In the usual LQC treatment one uses cylindrical functions, which are related with the connection variables through 
the holonomies (see e.g. 0[T3]). Here, we mimic this method for the spherical symmetric theory. First, we set up a 
polymeric Hilbert space "Hpoiy with a spin network like basis and define the action of the operators E v and K v on 
this state. We will then proceed by defining the gravitational part of our trial states centered around the classical 
Schwarzschild solution. 

In order to use the uniform discretization method, we consider a discretized setting. The symmetry reduced spatial 
manifold in the spherically symmetric case is a dimensional line (the radial direction). In analogy with LQC, we 
consider the bulk Hilbert space for the gravitational sector as 

Hpoly = L 2 (® N R Bo hr, ®Ndp ) , (26) 

where i?Bohr is the Bohr compactification of the real line, dpo is the Haar measure and N is the number of cells in 
the discretization. In this framework, the basis for a fixed graph g composed by N edges and vertices vj (with 
j = l,...,N) is: 



3 



(27) 
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with Hj e K. The variables satisfy the classical Poisson bracket 

{K v (i),E«>(j)} = G5 ij , (28) 
which suggests defining E v as a derivation operator 

^'-'isijry (29) 

Its action over the elements \p) of the basis of the gravitational Hilbert space is defined by 

E<o(i)\ft)=£ 2 p fi i \(t). (30) 

We associate the K v to point holonomies, and their action is defined by 

exp(ipK v (i))\fi) = \p + pei). (31) 

We want to remark that the Hilbert space we are using is a direct product of the Hilbert space of LQC for every 
lattice position i, which enables us to use techniques developed in LQC for our purpose. 
Finally, we construct the trial states as 



!•»> - tf 1 *> - n fk> «*> (-;bj ("< - ) • < 32 > 

with /Ltj centered at E v (i) — l 2 ^ = en(i), where the value of rx(i) will be determined in the following. 

Since the classical counterpart of the vacuum solution corresponds to vanishing scalar fields, we can ignore the matter 
part of the Hamiltonian H matt and focus on the gravitational part H vac . Demanding H mc = and additionally using 
the gauge K v = one gets 

-x(l-2A)+ ^ + g 3 V = => E-- 1 X + a 



(E^) 2 J ~ x/T^2A h a-C ' 

y x+a 

where C is an integration constant. In order to recover the Schwarzschild metric for A = 0, we set C = 0. Then we 
obtain 

r,(i)= 1 . (33) 

We find that, in contrast to [3], a is not the Schwarzschild radius Rs but rather 

a = R s Vl - 2A. 

As it was commented before, the introduction of the cosmological constant induces a rescaling of the radial coordinate. 

B. Matter part of the trial states 

In order to construct the matter part of the trial states we define an effective Hamiltonian as the expectation value 



of the matter Hamiltonian constraint over the gravitational part of the trial states (32). We compute the equations 
of motion corresponding to this effective Hamiltonian and get a differential equation for the matter variables <f> and 
. Working in the Fourier space, this equation turns out to be a Sturm-Liouville problem. We can construct the 
associated creation and annihilation operators without knowing the explicit solution to the equation (making use 
of the tools of the Sturm-Liouville theory). It is thus straightforward to introduce the Fock- vacuum state, which 
completes the construction of the trial state. 

Due to the usual factor ordering ambiguities in the quantization procedure, we need to choose one prescription for 
it. Here, we use the factor ordering of [jQ, namely putting the operators E v and P^ symmetrically around K v and <f> 
respectively. In this way, the discretized matter Hamilton constraint operator can be written as: 

v ' (E^ii)) 2 

r(i) r~~ . nr- 1 - 1 ( 34 ) 

- 2^VPHcj)(i + 1) - 4>(i)) VP*^= zm{ P K v {i))-= ~ ep vac . 

P y/Ef V Ev 
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Using the expectation values ( A3 1 and ( A9 ) we obtain the effective Hamiltonian constraint 

2 



H^tt = <**|&n» tt |**) - ( (P+) + r 4 {$') j f(r) - P - 



(35) 



where we went back to the continuum theory because it is easier to solve the corresponding differential equations than 
the original difference equations. 



The equations of motion for the effective Hamiltonian constraint ( 35 1 are: 



<f>(r,t) 



6P4> 

X ffeff 

P*(r,t) = jf^- 



2/(r)P*(r,t), 
= (2/(r)rY(r,t))' 



Now, we consider the Fourier like transformation 

<j>(r,u) 



0(r,t) 



'2u 



(e %ut C{w) + e- luJt C(uj)) du , 



P*(r t) - 1 ^M) 



(r, w) (e^^w) - e- <w *C(w)) . 



(36) 
(37) 

(38) 
(39) 



where the transformation for the momentum follows from (36). This leads to the following Sturm-Liouville differential 
equation 

(40) 



{B{r)(j>'(r, w))' + Lj 2 A(r)cj)(r, u) = , 



with 



4(r) = 



2/(r 



B(r) = 2/(r)r 4 



(41) 



According to the Sturm-Liouville theory, solutions of such differential equations form a basis of the Hilbert space 
L 2 (R, A(r)dr). This means that properly normalized functions <f>(r,uj) satisfy the orthogonality and closure relations 



dr A(r) <j>(r, u)4>(r, u) — 6 (lu — oj ) 
1 

du) 4>(r, uj)<p{r', u) = -jfyj S ( r - r ') ■ 
Consequently, the field and momentum operators are constructed as: 



(r,t 



(e lwt C' t (w) + e~^*C , (w)) du, 



P*(r,t)=iA(r) / ^^'^(^H^^))^- 



(42) 
(43) 

(44) 
(45) 

Considering the commutation relation C(u), C'(ui') = 5(uj — uj') and making use of the closure relation, we obtain 
the standard commutator for the field and its conjugate momentum operators 

(r,t),P (/,t)j = iS(r-r') . 

As it is inferred by the previous commutation relations, the operators C and C* act as annihilation and creation 
operators respectively. 

Finally, we construct the trial state as the direct product of the Gaussian state ( 32 ) and the Fock vacuum of the 
matter part: 

|v[4 rial ) = |^) ® |0) , (46) 

where |0) is the state annihilated by the annihilation operator C. In this treatment, by construction, the case of 
entanglement between matter and gravitational variables is excluded, because the split between the gravitational and 
the matter sector works only with the vacuum state. For excited states this requirement should be relaxed. 



IV. EXPECTATION VALUE OF THE MASTER CONSTRAINT 



In this section we use the trial state (46) to calculate the expectation value of the master constraint. First, wc 
construct the discrete master constraint out of ( 24 ) . We are able to write the principal order of the expectation value 
in orders of the lattice spacing. It shows that the master constraint is becoming small for bigger lattice spacings. 
Since the approximations made break down at a certain point one concludes that there is a minimum of the master 
constraint. However, the minimum value is not reached in the limit e — > 0. This fact suggests that there is no 
continuum limit. 



A. The discrete master constraint 



Using the uniform discretization technique we construct the master constraint ([2| as 

1 f H(r) 2 

H= - / dr^^i p , (47) 

2 J rEv p ' v ; 

where we used our gauge and took into account that the Hamiltonian is a density of weight one, so we added the 
determinant of the 3-metric in order to make contact with the continuum theory. The Planck length £ p ensures that 
H is dimensionless. Notice that the master constraint is proportional to 1/E V . This implies that the vacuum state 
corresponds to the zero loop state (p = 0) which is degenerate in the polymer representation. 

Nevertheless, the Hamiltonian constraint can be rescaled by a scalar function of the canonical variables with- 
out changing the first class nature of the constraint algebra, so we can solve this problem by rescaling H(x) — ?> 
H(x)^/E i P/(E x y. Then, in the discrete theory we obtain 

i i 

For computational purposes, we proceed to separate "matter" and "gravity" operators acting on their corresponding 
Hilbert spaces, i.e., we write the Hamiltonian constraint (24 1 as 



where 



H(i) = H£ln(i)ci(i) + ^nittW^) + H^ att (i)c 3 (i) + c 4 (z) , (49) 

H£L(i) = e p {eP^f + er^fm + 1) - 0«) 2 ) , (50) 

H [ 2tS) = il {-2r{i)pf(<P(i + 1) - 4>ii))) , (51) 

tfmittM^pPvac, (52) 



and 



*«) =- (l-2A)e- Ui + l)^4±m _ r(0 «"Vv<<->n ' U . I ,V ,,;,V 



(53) 

p 2 ' p 2 ) ' Ev>(i + 1) 2 E^(i) 2 ' 

In terms of these operators, we can write the expression for the discretized master constraint in the following form: 



+ 2c 12 (i)Hi]Ui)H£lt(i) + 2c 13 (i)H^ tt (i)H^Ui) + 2c 23 {i)H^l tt m { H tt (i) (54) 
+2cu(t)fT<&t(i) + 2 C24 «tf 1 ( n 2 i tt « + 2c 34 (i)H^ tt (t) + c 44 (i)) , 



where c jk (i) = Cj(i) ■ c k (i). 

As the matter and the gravitational sector are not entangled, i.e., (0|cjk(i)|0) = Cjk{i) and ff\H^l tt (i)\^/ 
flniattW) we can accomplish the computation of the expectation values separately. 
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B. Minimizing the master constraint 



In this part we treat the problem of minimizing the master constraint. All the intermediate calculations are included 
in the appendices. For the sake of simplicity, we go back to the continuum theory and assume that a is independent 
of the lattice position, which means: 



r(i) — > r , r(i + 1) — >• r + e . 
cr(i) — > a , a(i + 1) — > a . 



(55) 



Because of the independence of a with respect to the position we also set a = <7q . 

p 

As it was explained in [T], the approximation we considered in order to handle the expressions (neglecting higher 
powers of e) is inadequate for large values of e. More specifically, the approximation is valid up to values of e ~ 
10 -23 cm. In this range it is convenient to sort the terms in orders of 0((£ p /e) k e). Then, the non vanishing components 
of the master constraint ( 55 ) up to the principal order are: 



(cii(^itt) 2 ) = (3AM 4 / 2 + 3r 8 / 2 + 2A(ry (2/ 3 2 + hh)) 
(c 13 H£l tt H£l tt )=(A(r) 2 I* + r i I: 
(c 14 ^ tt ) = (A(r) 2 7 1 2 + r 4 7 2 ) 



-^A + offe 6 



1 1 



1 — exp 



<C2 2 (J?^ tt ) 2 ) 



= r 2 A(r) 2 (2/f + J 1 / 2 



1 

v 



exp 




<C33(/4?itt) 2 >=4A : 



\C-34-n matt 13 matt 



2 e 2 . 



1 



(C44) 



v 



2p 2 
1-2 exp I 



1 — exp 




a e 



exp 



(56) 

(57) 

(58) 

(59) 
(60) 
(61) 

(62) 



where Ji , I2 and ^3 are integrals coming from the computation of the expectation value of the Hamiltonian matter 
constraint and are given by 



h 



du>u)((j>(r, u>)Y 



1 



^-(0'(r,w)) 2 , 



la 



du (f)(r, uj)4>'{r, lo) 



(63) 



where (f>(r,u)) is the solution of the Sturm-Liouville equation (40). In the appendix [C] we give the expressions (C25) 
of these integrals for an approximate solution (f>o(r,uj) of the zeroth order of the Sturm-Liouville equation (B5). 
Finally, the main order of the expectation value of the integrand of the master constraint takes the form 



(HM)=^ T ifl-2exp^ V ^ 



e 16p 4 



cr e' 



exp 



(T e n 



O 



(64) 



Notice that the leading order term does not depend on the variable a. Also, since the equations are lengthy and they 
do not provide any further conclusions, we leave away higher order corrections. 

We are now in conditions to study the minimum of this master constraint. In figure [l] we plot the expectation value 
(for values cto = 10, n = 2 ) with respect to the lattice spacing e in the region of interest. We observe that the master 
constraint drops very fast for lattice spacings larger than the Planck scale. As in the case studied in jl], because of 
the break down of the approximation for e > 10~ 23 cm, we expect the master constraint to increase again. So, we can 
conclude that there is a minimum around e w 10~ 23 cm. 

We find that varying the variable a, which is connected to the Schwarzschild radius of the classical solution where 
the trial states are centered, has influence on the master constraint only at the subleading orders. So we can conclude 
that, as in the Minkowski case studied in pQ, we can construct the vacuum state for our minimally coupled system. 
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FIG. 1: Expectation value of the integrand of the master constraint for ctq = 10, n = 2 and p — 1 as a function of the lattice 
spacing e in the region near i p . 



In this paper we worked out an approximation for the vacuum state of a scalar field coupled to gravity with 
spherical symmetry. We focused on the Schwarzschild spacetime for the gravitational sector and we employed the 
method presented in pQ. 

The vacuum for the coupled system is given by the direct product of the Fock vacuum state for the scalar field 
and a Gaussian centered around the classical Schwarzschild solution for the gravitational sector. In order to deal 
with the dynamics, the uniform discretization technique is used. This setting allows us to develop a minimization of 
the (discrete) master constraint of the system using a variational method. In order to accomplish this we need to 
construct the Fock vacuum of the matter sector. We get this vacuum solving the Hamilton equations for an effective 
Hamiltonian, which turn out to be a Sturm-Liouville differential equation. It has to be pointed out that we did not 
manage to find an exact solution for this equation in the case we are interested in. Nevertheless, a suitable approximate 
solution is provided. Once the trial state is established, it is possible to perform the variational method in order to 
minimize the master constraint. Finally, we get an expression for the principal order of the expectation value of the 
master constraint with respect to the lattice spacing e and the Planck length t p . We found that, at least at the main 
order, the expectation value of the master constraint does not depend on the classical Schwarzschild radius. 

At this point it is worth comparing our results with the ones of Gambini, Pullin and Rastgoo. First, we have, at 
the leading order in the expectation value of the master constraint the same situation as in the cited paper, where 
the problem was worked out around a flat Minkowskian spacetime. As in that case, we also have a diverging master 
constraint in the continuum limit, but the theory gives a good approximation to general relativity for small values of 
the lattice separation. The reason for that, clearly explained in [I], is due to the gauge fixing procedure employed to 
get rid of the diffeormorphism constraint at the classical level. 

This work can be extended in several ways. For example, performing a polymeric quantization for the coupled sys- 
tem, or trying to work out this computation avoiding the gauge fixing of the diffcomorphism constraint. Nevertheless, 
these options would clearly rise the technical complexity of the problem. Anyway, in our opinion, this is a promising 
research line. A deeper understanding of this setting will provide us with a suitable vacuum state in this context 
opening the possibility of performing a detailed study of the Hawking radiation within the LQG framework. We left 
this problem for future investigation. 



V. 



CONCLUSIONS 
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Appendix A: Computation of expectation values for basic operators 

In this appendix we give the expressions for several basic expectation values that are used along the paper. 
First, we compute the expectation value of the operator E v (i): 



where \^s) is the gravitational part of the trial state given by the equation (32 1 



For the expectation value of trigonometric functions we proceed in the same way, and find 



;*, T |cos(2 /) AV(/))|v^ T ) = ( -J' ) . ,A1.) 



,2 



(^|cos(4pif v (i))|^) = exp^--^J , (A2) 
(*g\— f ^=ffa{j>k v (%))— F L=\* 3 ) = 0, (A3) 



(tta U sin(p^(z)) > = 0, ( A 4 ) 

(E^(t)) 3 / 2 (Ef(i)) 3 / 2 

,*>g\ * sin(3p£ y (t))— = 0. (A5) 

'Ev(i) yEf{i) 



Now we use the inverse (E v ) 3 / 2 operator, following the prescription given for LQC in [TSJIIS]. We obtain for the 
eigenvalue of \fl): 

(E^)y^\fl) = ^ + p) 3/4 - Ok - p) 3/4 ) 6 im) . (A6) 

This can be used to calculate expectation values of operators involving the inverse of 

(\ 2m/3 
^rrr- I |**>, (A7) 

where m > k > can be chosen arbitrarily, depending on the prescription taken for applying the Thiemann's trick. 
In fact, there is the same kind of ambiguity in LQC for the expression of inverse volume operators. 

We now concentrate on lattice spacings re > !p. In this regime we obtain the following approximation 
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As commented before, the dependence of ( A8 ) on m shows an ambiguity which occurs due to the use of the inverse 
operators. More specifically, for the inverse operators needed in our case we obtain: 



1 



1 



1 



e 2 ri (i)2 eVi(i) 4 V > V 



(£*>(»))* e 4 ^) 4 ^(ir^ p ^ 6 



(A9) 
(A10) 



with 



a(z) 



( 5m a 2,3 

O H (Til) 

I 24 4 W 



where TO a > 2 and > 4. Analogously one gets 



**|^— cos(2pK ip (i))-^— \V 9 ) =exp ( 



E*{i) 



a{i) 



e 2 ri(z) 2 e 4 ri(z) 4 ' w p / I e 6 



(All) 



with 



7(«) 



5m 



12 



1 P 2 + 



m 7 > 1. 



Appendix B: Solution to the Sturm-Liouville problem 



We consider the Sturm-Liouville problem (40) for the leading order <f>o(r,ui) in the expansion of the solution in 
powers of life 2 . 



(2L (l - ~) rVoM)' + ^^Mr,") = 0, 

r 

where L = (1 — 2A). We were able to find the function 



(Bl) 



(B2) 



where 



a loa 



- 1 



is the usual tortoise coordinate of the external region of the Schwarzschild metric. The approximate solution (B2) 
fulfills 



(2L (l - ~) r 2 ^(r, W ))' + ^^^(r.u) = -2^L0 o (r, W ) , 



(B3) 



which for r 3> a properly approximates the Sturm-Liouville equation (Bl|. We can observe in figure [2] that indeed 
the approximation is very accurate (it overlaps completely the exact numerical solution). Another possible option 
could be motivated by the fact that in the regime r»a the equation (Bl) becomes the Sturm-Liouville problem of 
the Minkowski case, studied in [I], with exact solution 



i(r,w) = - sin 



( w 
\2L 



(B4) 



However this Minkowskian solution does not succeed in approximating the exact n umer ical solution, as it is illustrated 
by figurepl We will therefore work with (f>a(r, ui). As mentioned before solutions of (Bl ) should fulfill the orthogonality 
relation (|42[). We can use this fact to normalize our solution, which is then given by 



2 1 



\2L 



(B5) 
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FIG. 2: Comparison of the approximate solution (f>o(r,uj) (solid line) given by equation (B2l with the numerical one (dashed 
line) for uj = L/a. Notice that they overlap completely, showing the quality of the approximation. 




FIG. 3: Comparison of the Minkowskian solution <f> m (r, u) (solid line) given by the equation (B4l with the numerical one 
(dashed line) for uj = L/a. Notice that 4>m{r, w) is not a suitable approximation. 



Appendix C: Coefficient operators and expectation values 
1. Coefficient operators 



The operator form of the coefficients of the discrete master constraint (24) is: 
1 

Cll 



Cl2 



Cl3 



C22 



C23 



1 sin(pK v (i)) 1 
(£*>(i))3/2 p (£V(;))3/2 ' 

1 



J- Ui 

V l E*{i) 2 Ef(i) 



cos(2pK v (i)) 



Ef{i) 



1 sm(pK v (i)) 1 
' {E*{i)yl 2 P (^(i))V2 



(CI) 
(C2) 
(C3) 

(C4) 

(C5) 



14 



C33 = 1 



(C6) 



cm = \ r{i)—-r{i + l) 



I - cos(2 P k v {i + I)) 
2p 2 



r{i) 



1 cos(2pK^(i)) 1 



- (1 - 2A)e 
r^ + l) 3 e 2 



1 

r(i) 3 e 2 



(C7) 



E<P(i) 2p 2 Ev(i) E<P(i + l) 2 Ev>(i) 2 E<P(i)± ' 



C24 



-(i + 1) 



3,2 



-r(i + l) 



1 - cos{2pK v (i + 1)) 



r(i) 



£^(i + l) 2 2p 2 

1 3sin(p.fir<p(i)) - sin(3pJ¥' ¥ ,(i)) 1 



(l-2A)e 



E*{iyi 2 



4p 3 



Ev(iy/ 2 



-r(i) 



1 sin(pK'y(i)) 1 

(l^jjv 2 p {Ef{i)y/ 2 

3 2 1 &n{pK<p(i)) 1 



(£¥>(j))3/2 



(^(i))3/2 



C34 = r(i + 1) 



1 - cos(2pKJi + 1)) 1 - cos(2pKJi)) r(i + l) 3 e 2 r(i) 



\3 C 2 



2p 2 



r(i) 



2p 2 



Ev(i + 1) 2 Ef(iy 



+ (l-2A)e, 



(C8) 



(C9) 



c 44 = (1 - 2A) 2 e 2 + 2(1 - 2A)e r(i + 1) 



1 - cos(2pK ip (i + 1)) 1 - cos(2pK ip (i)) 



2p 2 



-r(i)- 



2p 2 



n2 3 - Acos(2pK ip {i + 1)) + cos(4pK v (i + 1)) 
^p^ 

+ r (i) 2 lJ ~ ^^^p^v)) + cos(4pify w/ + 2 , ^ 



j 3 - 4cos(2p^(i)) + cos(4pi^(i)) 



-2r(z)r(i + l) 2 ^ 2 



1 - cos(2pif v (i)) r(i+l) 



2p 2 



2p 2 



(1 -2A)e 



2p 2 
r{i + l) 3 e 2 
Ev(i + l) 2 



+ 2 \r(i + l) 



1 - cos(2 P K v (i + I)) _ Ki) \ r(i) 3 e 2 4 2 1 cos(2pir y (i)) 1 

2p 2 2p 2 1 ' J E*(i) 2 1 J Ef{i) 2p 2 Ef{i) 



4,2 



2r(i + l) 4 e 



1 cos(2piL,(i + l)) 1 r(i + l) 6 e 4 _r(i) 3 e 2 r(i + l) 3 e 2 r(i) 6 e 4 



+ l) 



2p 2 



Ev(i + 1) £V(i + l) 4 Ev(i) 2 £¥>(£ + 1) 2 £:v(j)4 



(CIO) 



Using (A1)-(A5) and (A9)-(A11| one can now calculate the expectation values of the coefficients (c-ij) 
(*? ial |ciiT^f ia 7 = {^ff\cij^9hWe remark that 



(ci 2 ) = (C23) = (qm) = . 



(Gil) 



2. Expectation values of the "matter Hamiltonians" 

Let us focus now on the expectation valu es of the expressions of the discrete master constraint which contain 
"matter Hamiltonians". Since the coefficients (Cll) vanish, we only need to compute ((H^l tt (i)) 2 } , ((H^l tt (i)) 2 } and 
(H^l tt (i)). The continuum version of (501 can be written as 



H {1) -e 2 



£ 2 (pV,*) 2 +''VM) 2 ) • (C12) 



Using {H^l u (i)) = (0\H^ ] att {i)\0) , (0\C(u})C^ (uj')\0) = 6(u - w') and that all the other combinations of the C- 
operators vanish, one gets: 



r(l) \ _ 6 *P ( ^„\2 / j. . .^2 , „4 / j. . 1 /I/ 



(^tt) = Yr (r) X d ""(Mr,u)y+r* du-{Ur,u)Y) ■ (013) 
For (i?^ tt ) 2 , with the given factor ordering we obtain 

(i?W tt ) 2 = A 4 (i*(r, tf + 2r 4 P^, t)j>'(r, t) 2 P+(r, t) + r 8 0'(r, i) 4 ) . (C14) 
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The expectation value of (C14) becomes 



where 



((^2tt) 2 ) = ^ (ZMrfll + 2A(r)V (2/ 8 + hh) + 3r 8 / 2 2 ) , (C15) 



/, = / '^1. : . n (,,ir. /, = / ^(^(r,^)) 2 73=/ I-',!,! r.^i. (CICi) 



Now, for (-ff^att) 2 ! which following from (51) has the continuum limit 

(H^luf = rXe 2 PHr,")4>'(r,u;fP*(r,u>) , (C17) 

we obtain 

((i^ tt ) 2 ) = r 2 A(r) 2 , 2 l% (2/ 3 + hh) . (C18) 

Summarizing and going back to the discrete theory we get the following expectation values for the "matter Hamilto- 
nians" : 

(A£Ui)) = f ^ {A{r{i)fh{i) + r(i)%(i)) , (C19) 

((A ( nittW) 2 > = | e 6 (3^(^W) 4 /iW 2 +3r(z) 8 / 2 ( i ) 2 + 2A(r(z))r( i ) 4 (2/ 3 ( l ) 2 + / 1 W/ 2 (z))) , (C20) 

((^i 2 2 tt (i)) 2 )=r(i) 2 A(r(i)) 2 e 4 ^(27 3 (i) 2 +/ 1 (i)/ 2 (i)) , (C21) 

(ff£2 tt (i))=2Ae, (C22) 

((^ tt (i)) 2 )=4A 2 e 2 , (C23) 

(H£l«H { *lM = Aiy (A{r{i)fh{i) + r(i) 4 7 2 (i)) , (C24) 

where we used p vac = 2A/£ 2 . 

Now all we need to calculate to get all the components of the master constraint (55) is the three integrals (C16). 
Using the approximate solution (B5) and taking into account the discretization of the integral one gets: 

2 1 /tt 2 L 2 cosf^fr*) , an(Hr") \ , s 



Trr 2 \e 2 2(r*) 2 2(r*) 2 



I I ( „.{2ir „\ . , {2* 2 1 1 sin 2 (^r* 



h = —Al-Gi[—r*) +log — r 



7r r 4 \ \eL J \eL J J tt r 3 1 

-L 2 + ^ 



1 1 ^ 2 L 2 cos(ff r*) - sin(f£r*) 



2ttL 2 ^2 (1 _ ay \e 2 2(r*) 2 2(r*) 2 ~ < /■' 



(C26) 



11 1 /sin(^r*) 2 cos( 2 fr*) \ 2 1 /tt sin^r*) , 
/3 = ^L^1^|( V 2(r*) 2 ^ " cr* ~ *^7 ~ 2r* L) ' ' C -' 1 

where Ci(a;) is the cosine integral function 2 and 7 is the Euler's constant 
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